A Renormalization Group Approach for Highly Anisotropic Fermion systems 
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I apply a two-step density-matrix renormalization group method to the anisotropic two- 
dimensional tight-binding model. This study, which is a prelude to the study of models of quasi-one 
dimensional materials, shows the potential power of this approach for anisotropic fermionic models. 
I find a ground-state energy which agrees with the exact value up to four digits for systems as large 
as 24 X 25. This open new opportunities for simulations of fermions in two dimensions. 



Quasi-one dimensional organic and inorganic 
materials have been the object of an important theo- 
retical interest for the last three decades. The essential 
features of their phase diagram may be captured by the 
simple anisotropic Hubbard model (AHM), 



H = Y.^4l^aCi+l,l,<r + h.C.) 

-t±_ ^{c\i^„CiJ+i^cr + h.C.) + 

id id, a 
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or a more general Hubbard-like model including longer 
range Coulomb interactions. For these highly anisotropic 
materials, t± ^ Over the years, the AHM has re- 
mained a formidable challenge to condensed-matter the- 
orists. Some important insights on this model or its low 
energy version, the g-ology model, have been obtained 
through the work of Bourbonnais and Caron 0, Q and 
others. They used a perturbative renormalization group 
approach to analyze the crossover from ID to 2D at low 
temperatures. More recently, Biermann et al. 5] applied 
the chain dynamical mean-field approach to study the 
crossover from Luttinger liquid to Fermi liquid in this 
model. Despite this important progress, crucial informa- 
tion such as the ground-state phase diagram, or most no- 
tably, whether the AHM displays superconductivity, are 
still unknown. So far it has remained beyond the reach 
of numerical methods such as the exact diagonalization 
(ED) or the quantum Monte Carlo (QMC) methods. ED 
cannot exceed lattices of about 4x5. It is likely to re- 
main so for many years unless there is a breakthrough 
in quantum computations. The QMC method is plagued 
by the minus sign problem and will not be helpful at low 
temperatures. The small value of t±_ implies that, in or- 
der to see the 2D behavior, it will be necessary to reach 
lower temperatures than those usually studied for the 
isotropic 2D Hubbard model. Hence, even in the absence 
of the minus sign problem, in order to work in this low 
temperature regime, the QMC algorithm requires special 
stabilization schemes which lead to prohibitive cpu time 



I have shown in Ref.0 that this class of anisotropic 
models may be studied using a two-step density-matrix 
renormalization group (TSDMRG) method. The TS- 
DMRG method is a perturbative approach in which the 
standard ID DMRG is applied twice. In the first step, 
the usual ID DMRG method is applied to find a set 
of low lying eigenvalues e„ and eigenfunctions |(/)„) of a 
single chain. In the second step, the 2D Hamiltonian is 
then projected onto the basis constructed from the tensor 
product of the |(/)„)'s. This projection yields an effective 
one-dimensional Hamiltonian for the 2D lattice. 



H 
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where Eni^] is the sum of eigenvalues of the different 



chains, En 



|<f>ll[„]) are the corresponding 



eigenstates, |*||[„]> = |^„i)|^„2)---l'/'nz,); c],p and 
the renormalized matrix elements in the single chain ba- 
sis. They are given by 



(3) 
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where n.i represents the total number of fermions from 
sites 1 to z — 1. For each chain, operators for all the sites 
are stored in a single matrix 
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Since the in-chain degrees of freedom have been inte- 
grated out, the interchain couplings are between the 
block matrix operators in Eq. ( ^ which depend only 
on the chain index I. In this matrix notation, the effec- 
tive Hamiltonian is one-dimensional and it is also studied 
by the DMRG method. The only difference compare to 
a normal ID situation is that the local operators are now 
77152 X TTIS2 matriccs, where 77752 is the number of states 
kept during the second step. 

The two-step method has previously been applied 
to anisotropic two-dimensional Heisenberg models. In 



Ref. [Sj , it was applied to the t — J model but due to the 
absence of an exact result in certain limits, it was tested 
against ED results on small ladders only. A systematic 
analysis of its performance on a fermionic model on 2D 
lattices of various size has not been done. In this letter, 
as a prelude to the study of the AHM, I will apply the 
TSDMRG to the anisotropic tight-binding model on a 
2D lattice, i.e., The AHM with U = 0. 1 perform a com- 
parison with the exact result of the tight-binding model. 
I was able to obtain agreement for the ground-state en- 
ergies on the order of 10""* for lattices of up to 24 x 25. 
I then discuss how these calculations may extend to the 
interacting case. 

The tight-binding Hamiltonian is diagonal in the mo- 
mentum space, the single particle energies are, 
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FIG. 1: Low-lying states of the ID tight-binding model (full 
line) and of the ID Heisenberg spin chain (dotted line) for 
L — 16 and ms2 = 96. 



efe = —2t\\coskx — 2t±cosky — ^, (7) 

with k = {kx, ky), kx — UxTt / (Lx + l) and ky — nyTT/{Ly+ 
1) for open boundary conditions (OBC); Lx, Ly are re- 
spectively the linear dimensions of the lattice in the par- 
allel and transverse directions.. The ground-state energy 
of an N electron system is obtained by filling the low- 
est states up to the Fermi level, E[q]{N) — X]fc<fcp ^fe- 
However in real space, this problem is not trivial and it 
constitutes, for any real space method such as the TS- 
DMRG, a test having the same level of difficulty as the 
case with U 0. This is because the term involving U is 
diagonal in real space and the challenge of diagonalizing 
the AHM arises from the hopping term. 

I will study the tight-binding model at quarter filling, 
N/LxLy — 1/2, the nominal density of the organic con- 
ductors known as the Bechgaard salts. Systems of up to 
Lx ^ Ly L X 

(L+1) = 24x25 will be studied. During the 
first step, I keep enough states (msi is a few hundred) so 
that the truncation error pi is less than 10^^. I target the 
lowest state in each charge-spin sectors Afa;±2, Nx±l, Nx 
and 5*2 ± 1, 5*2 ± 2, Nx is the number of electrons within 
the chain. It is fixed such that Nx/Lx = 1/2. There is a 
total of 22 charge-spin states targeted at each iteration. 

For the tight-binding model, the chains remain discon- 
nected if tj_ < eo{Nx + l) - eo{Nx) or tj_ < eo{Nx) - 
£oiNx ~ 1), where iV^; is the number of electons on sin- 
gle chain. In order to observe transverse motion, it is 
necessary that at least t± > eo{Nx -f 1) — eo{Nx) and 
t± ^ eo(-^a;) — eo{Nx — 1)- These two conditions are sat- 
isfied only if fi is appropietly chosen. The values listed in 
Table (jH) corresponds to fi = {eo{Nx + l)-€o{Nx-l))/2. 
This treshold varies with L. I give in Table iQJ the val- 
ues of t± chosen for different lattice sizes. In principle, 
for the TSDMRG to be accurate, it is necessary that 
Ae = e„^ — eo, where e„^ is the cut-off, be such that 
Ae/t± ^ 1. But in practice, I find that I can achieve 
accuracy up to the fourth digit even if Ae/t± « 5 using 
the finite system method. Five sweeps were necessary 





8x9 


16 X 17 


24 X 25 


t± 


-0.28 


-0.15 


-0.1 




-1.2660 


-1.3411 


-1.3657 


Ae/tx 


6.42 


5.40 


5.78 



TABLE I: Transverse hopping and chemical potential used in 
the simulations for different lattice sizes 

to reach convergence. Note that this conclusion is some- 
what different from my earlier estimate of Ae/tj_ « 10 
for spin systems 01 • This is because in Ref.(0]), I used 
the infinite system method during the second step. 

The ultimate success of the TSDMRG depends on the 
density of the low-lying states in the ID model. For 
fixed ms2 and L, it is, for instance, easier to reach 
larger Ae/Jj_ in the anisotropic spin one-half Heisen- 
berg model, studied in Ref.0, than Ae/t± for the tight- 
binding model as shown in Fig^ For L = 16, ms2 — 96, 
and J^^ = tjL = 0.15, I find that Ae/J±_ « 10, while 
Ae/ij. « 5. Hence, the TSDMRG method will be more 
accurate for a spin model than for the tight-binding 
model. Using the infinite system method during the 
second step on the anisotropic Heisenberg model with 
J± = 0.1, I can now reach an agreement of about 10^^ 
with the stochastic QMC method. 

Two possible sources of error can contribute to reduce 
the accuracy in the TSDMRG with respect to the conven- 
tional DMRG. They are the truncation of the superblock 
from 4xTOSi states to only ms2 states and the use of three 
blocks instead of four during the second step. In Table 
(|IH) I analyze the impact of the reduction of the number of 
states to ms2 for three- leg ladders. The choice of three- 
leg ladders is motivated by the fact that at this point, 
the TSDMRG is equivalent to the exact diagonalization 
of three reduced superblocks. It can be seen that as far as 
t± > eo{Nx+l)-€o{Nx) andii > eo(iV,)-eo(iV.-l), the 
TSDMRG at this point is as accurate as the ID DMRG. 
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ms2 8x3 16 x3 24 x3 



64 


-0.241524 


-0.211929 


0.204040 


Exact 


-0.241524 


-0.211931 


0.204049 


TABLE II: 


Ground-state 


energies of three-leg ladders. 


ms2 


8x9 


16 X 17 


24 X 25 


64 


-0.24761 


-0.21401 


0.20504 


100 


-0.24819 


-0.21414 


0.20509 


120 


-0.24832 


-0.21419 




Exact 


-0.24857 


-0.21432 


0.20519 



TABLE III: Ground-state energies for different lattice sizes; 
a single state was targeted in the second step. 

Note that the accuracy remains nearly the same irrespec- 
tive of L as far as the ratio Ae/t± remains nearly con- 
stant. Since Ae decreases when L increases, t± must be 
decreased in order to keep the same level of accuracy for 
fixed ms2- In principle, following this prescription, much 
larger systems may be studied. Ae/t± does not have to 
be very large, in this case it is about 5, to obtain very 
good agreement with the exact result. 

The second source of error is related to the fact that 
the effective single site during the second step is now a 
chain having ms2 states, I am thus forced to use three 
blocks instead of four to reduce the computational bur- 
den. In Table Ijllip . it can be seen that this results in a 
reduction in accuracy of about two orders of magnitude 
with respect to those of three leg-ladders. These results 
are nevertheless very good given the relatively modest 
computer power involved. All calculations were done on 
a workstation. 

The DMRG is less accurate when three blocks are used 
instead of four. This can be understood by applying the 
following view on the formation of the reduced density 
matrix. The construction of the reduced density matrix 
may be regarded as a linear mapping w^r : F* —> E, 
where E is the system, F is the environment and, F* 
is the dual space of F. Using the decomposition of the 
superblock wave function \I>[o] = J2i 4'i ® 4>f'i with 0f G 
E and 4>f e F, for any (/)* eF*, 

u^m^Y.^'^icb^'Hf. (8) 

i=l 

Let \k), k = l,...dimE and \l), I = l,...dimF be the 
basis of E and F respectively. Then, |/) has a dual ba- 
sis {l*\ such that {l*\l) — Si J-'. The matrix elements 
of It* in this basis are just the coordinates of the su- 
perblock wave funtion $[0]^.^. The rank r of this map- 
ping, which is also the rank of the reduced density matrix 
is always smaller or equal to the smallest rank of E and F, 
r < Min{dim'E, dimF). Hence, if rns2 states are kept in 



ms2 


8x9 


16 x 17 


64 


-0.24803 


-0.21401 


100 


-0.24828 


-0.21417 


Exact 


-0.24857 


-0.21432 



TABLE IV: Ground-state energies for different lattice sizes; 
three states were targeted in the second step: the ground state 
itself and the lowest states of = and Sz ~ 1 sectors. 

the two external blocks, the number of non-zero eigenval- 
ues of p cannot be larger than ms2- Consequently, some 
states which have non-zero eigenvalues in the normal four 
block configuration will be missing. A possible cure to 
this problem is to target additional low-lying states above 
\I'[Q](iV). The weight of these states in p must be small 
so that their role is simply to add the missing states not 
to be described accurately themselves. A larger weight 
on these additional states would lead to the reduction 
of the accuracy for a fixed ms2- In table (|IV|) . I show 
the improved energies when, besides the ground state, 
I target the lowest states of the spin sectors Sz = —I 
and Sz = +1 with TV electrons. The weights were re- 
spectively 0.995, 0.0025, and 0.0025 for the three states. 
This lowers E[q] (TV) in all cases, but the gain does not 
appear to be spectacular. But I do not know whether 
this is due to my choice of perturbation of p or whether 
even the algorithm with four blocks would not yield bet- 
ter E[Q]{N). If the lowest sectors with iV -I- 1 and — 1 
electrons which have Sz = ±0.5 are projected instead, I 
find that the results are similar to those with Sz — 
sectors, there are possibly many ways to add the missing 
states. A more systematic approach to this problem has 
recently been suggested It is based on using a local 
perturbation to build a correction to the density matrix 
from the site at the edge of the system. Here, such a per- 
turbation would be Ap = acjpci, where a is a constant, 
awlO"'^ — 10~^, and cj, ci are the creation and annihi- 
lation operators of the chain at the edge of the system. 
This type of perturbation resulted in an accuracy gain of 
more than an order of magnitude in the case of a spin 
chain The three block method was found to be on 
par with the four block method. It will be interesting to 
see in a future study how this type of local perturbation 
performs within the TSDMRG. 

To conclude, as a first step to the investigation of inter- 
acting electron models, I have shown that the TSDMRG 
can successfuly be applied to the tight-binding model. 
The agreement with the exact result is very good and 
can be improved since the computational power involved 
in this study was modest. The extension to the AHM 
with U is straightforward. There is no additional 
change in the algorithm since the term involving U is lo- 
cal and thus treated during the ID part of the TSDMRG. 
The role of U is to reduce Ae as shown in Fig|21 For fixed 
L and ms2, Ae decreases linearly with increasing U. For 
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FIG. 2: Width Ae for the low-lying states of the ID Hubbard 
chain as function of Z7 for L = 16 and ms2 = 128. 



L = 16 and ms2 = 128, I anticipate that for C/ ^3 the 
interacting system results will be on the same level or 
better than those with ms2 = 100 for the same value of 
L. 
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